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Cd Abstract 

^ G^ We consider Ising models in d = 2 and d = 3 dimensions with 

nearest neighbor ferromagnetic and long-range antiferromagnetic in- 
^ teractions, the latter decaying as (distance)"^, p > 2d, at large dis- 

.^ tances. If the strength J of the ferromagnetic interaction is larger 

■^ than a critical value Jc, then the ground state is homogeneous. It 

^\^ has been conjectured that when J is smaller than but close to Jc the 

ground state is periodic and striped, with stripes of constant width 
^t h = h[J), and /i — )■ oo as J — t- J^ . (In d = 3 stripes mean slabs, not 

(T^ columns.) Here we rigorously prove that, if we normalize the energy in 

'^ such a way that the energy of the homogeneous state is zero, then the 

K*" ratio eo(J)/es(J) tends to 1 as J — )■ J~ , with es(J) being the energy 

k>< per site of the optimal periodic striped/slabbed state and eo(J) the 

^ actual ground state energy per site of the system. Our proof comes 

with explicit bounds on the difference eo{J) — es(J) at small but fi- 
nite Jc — J, and also shows that in this parameter range the ground 
state is striped/slabbed in a certain sense: namely, if one looks at a 
randomly chosen window, of suitable size i (very large compared to 
the optimal stripe size h{J)), one finds a striped/slabbed state with 
high probability. 
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1 Introduction and main results 

We consider Ising models in two and three dimensions on the square lattice 
with the formal Hamiltonian 

ff = -iEK-.-i) + E%^ (1.1) 

(x,y> {x,y} 

where the first sum ranges over nearest neighbor pairs in Z'^, d = 2, 3, the 
second over pairs of distinct sites in Z'^, and the exponent p is chosen to 
satisfy p > 2d, for reasons that will become clear below. For more general 
values of p, this model is used to describe the effects of frustration induced 
in thin magnetic films by the presence of dipolar interactions [p = 3) or in 
two-dimensional charged systems by the presence of an unscreened Coulomb 
interaction {p = 1) [D El gl El il El El El [IHl [111 Ell |22l [23l El ESI ESI E^ 
see also [3l[l2l[T3l[Tll[T5l[l6l[I71[l8]fora more detailed introduction to the 
subject, as well as for previous rigorous results. The competition between 
short range ferromagnetic and long-range antiferromagnetic interaction is 
believed to be responsible for the emergence of non-trivial "mesoscopic pat- 
terns" in the ground and low-temperature states of the system. Let us be 
more specific. As proved in [17], if J > Jc, with 

"^'■^ 5Z (^2 + 1 ±|2W2 ' (1-2) 

then there are exactly two ground states, (Tx = +1 Vx G Z,'^, and a^ = —1 
Vx G Z*^. Note that Jc is the value of the ferromagnetic coupling such that the 
energy of a straight domain wall configuration, i.e., a configuration consisting 
of half the spins minus (those at the left of a vertical straight plane) and half 
the spins plus (those at the right of the same plane), vanishes. If J ^ Jc, 
the ground state is certainly non-homogeneous. There is evidence that the 
transition to the ferromagnetic phase as J — )■ J^ takes place via a series of 
"microemulsion phases" characterized by phase separation on a mesoscopic 
scale that is large compared to the lattice and small compared to the scale of 
the whole sample; see e.g. [201 EHl ESI EH] for a discussion of this phenomenon 
in the case of Coulomb [p = 1) and dipolar {p = 3) interactions. More 
precisely, at zero temperature, the transition to the ferromagnetic state is 
expected to take place via a sequence of transitions between periodic striped 
or slabbed states, depending on dimensionality, consisting of stripes/slabs 
(either vertical or horizontal) all of constant width h{J) and of alternating 
sign. If we denote by es{h) the energy per site in the thermodynamic limit 
of periodic striped/slabbed configurations consisting of stripes/slabs all of 



size /i, the optimal stripe/slabs width h{J) can be obtained by minimizing 
es(/i) over /i G N, and turns out to be of the order (Jc — J)~p-d-i. Let us 
denote by es(J) := es{h{J)) the optimal striped/slabbed energy per site and 
by Co (J) the actual ground state energy per site in the thermodynamic limit. 
Our main result can be summarized in the following theorem: 

Theorem 1. As to J ^ Jc from below, we have 

lim ^ = 1 . (1.3) 

j^j- es(J) 



Eq. (1.3) is a strong indication of the conjectured periodic striped/slabbed 



structure of the ground state. The proof of Eq. (1.3) comes with explicit 



bounds on the speed of convergence to the limit, namely 

44 = 1 + o((j, - j)(.-i)r."-i)) . (1.4) 

It also comes with explicit bounds on the energy cost of the "corners" . This 
notion was introduced in ^17j for the two-dimensional case; every time that 
a domain wall bends by 90°, hence creating a corner (or an edge corner, as 
we call it, in three-dimensions: this is an edge where two plaquettes come 
together at 90°), we pay a positive energy cost, at least in the case that the 
corner density is sufficiently high. Combining this remark with our a priori 
bounds on the ground state energy, we find that the ground state has a 
density of corners that is smaller than (Jc — JY^^'^~^^: therefore, if we look in 
a random window of proper side C.' (much larger than the optimal stripe/slab 
width h{J) ~ (Jc — J)^p-'^~^ , and much smaller than the typical separation 
between corners ~ {Jc — J)~^^^'^'^^), the ground state restricted to such a 
window is striped/slabbed, with stripes/slabs of width close to the optimal 
size h{J). Our proof presumably adapts to any dimension, e.g., d = 10, 11 
or 26, and the interested reader can extend the arguments in Appendix |D] if 
desired. 

The logic of the proof goes as follows. We first derive an alternative 
representation of the energy in terms of droplet self-energies and droplet- 
droplet interactions. Next, for the purpose of a lower bound, we localize the 
energy into squares/cubes of side i (to be optimized over), and we show that 
the localized self-energy of every droplet with at least one corner along its 
boundary is positive; therefore, we can eliminate all such droplets, after which 
we are left only with striped/slabbed droplets. Finally, refiection positivity 
shows that the optimal striped/slabbed configuration is periodic. 



2 Droplets and self-energies 

Defining r := 2( J — Jc), the optimal periodic striped energy per site has the 
form: 

es(J) = -Cs(r)|r|(P-'^)/(P-^-i), (2.1) 

with Cs{t) = C3(0) + O(|r|2/(p-'^-i)) asymptotically for r -^ Q- , for a suitable 
Cs(0) > 0. This result follows from the explicit minimization of es{h), see 
Appendix |A| and can also be understood in terms of a balance between "line" 
or "plane" energies and line-line or plane-plane interactions, see [171 Section 
II]. We note that the computation in Appendix lA] also shows that the optimal 
stripe/slab width is 

h* = argmine,(/i) = C's(r)|r|'^/(P~^-^) , (2.2) 

with C's(r) = C'3(0) + O(|r|2/(p-"'-i)) asymptotically for r -^ 0", for a suitable 
6*3(0) > 0. Of course, eo(J) < es(J). Our purpose is to get a comparable 
lower bound, of the form 

eo{J) > -C3(0)|r|(P-'^)/(^-'^-i) (1 + 0{\Tf)) , (2.3) 

for some positive /3. The strategy borrows some ideas from those in pTl 
Appendix A]. 

From now on, for the purpose of simplicity of exposition, we restrict 
ourselves to two dimensions. We shall explain how to adapt the proof to 
three dimensions in Appendix [Dj We need to recall the definitions of contours 
and droplets. Let us first define the finite volume Hamiltonian for our system: 

^a(^a) = -JJ2 (^-^y - 1) + E %^^ + ^a(^aI^*) • (2.4) 
(x,y); {x,y}: 

x,yeA x,yeA 

Here A C Z^ is a square box, a^ = {axjxeA £ {il}^ is the spin configuration 
in A, a* = {o"*}xez2 G {±1}^ is a boundary condition and 

xeA, yeA=: xeA, yeA': ' "^ ' 

|x-y|=l 

In the discussion below, we shall consider + boundary conditions: this means 
that a* = {+1}^'. 

Given a^, we define A to be the set of sites at which a^ = —1, i.e., 
A = {xgA : o"x = — !}• Around each x G A we draw the 4 sides of the 



unit square centered at x and suppress the sides that occur twice: in this way 
we obtain a closed polygon r(A) which can be thought of as the boundary of 
A. Each side of r(A) separates a point x G A from a point y ^ A. At every 
vertex of r(A) fl (Z^)*, with (Z^)* the dual lattice of Z^, there can be either 
2 or 4 sides meeting. In the case of 4 sides, we deform the polygon slightly 
by "chopping off" the edge from the squares containing a — spin. See Figure 

m 



(a) 



+ 



(b) 



Figure 1: In the case that 4 sides of the closed polygon r(A) meet at a 
vertex v, we slightly deform r(A) so that the two squares containing a — 
spin become disconnected from the vertex itself. Case (a) represents the 
situation where the minus spins are located at NE and SW of v, before and 
after the "chopping" . Case (b) represents the situation where the minus spins 
are located at NW and SE of v, before and after the "chopping". 



When this is done r(A) splits into disconnected polygons Fi, . . . , F,. which 
are called contours. Note that, because of the choice of + boundary con- 
ditions, all the contours are closed. The definition of contours naturally 
induces a notion of connectedness for the spins in A: given x, y G A we 
shall say that x and y are connected if and only if there exists a sequence 
(x = xo, xi, . . . , x„ = y) such that x^,, x^+i, m = 0, . . . , n — 1, are nearest 
neighbors and none of the bonds (xm,Xm+i) crosses F(A). The maximal 
connected components 6i of A will be called droplets and the set of droplets 
of A will be denoted by T>{A) = {6i, ... ,6s}. Note that the boundaries T{6i) 
of the droplets Si G ^(A) are all distinct subsets of F(A) with the property: 

uLiF(50 = r(A). 

Given the definitions above, let us rewrite the energy H\{ajy) of a^ with 
+ boundary conditions as 



rer(A) SeviA) (5,5') 



(2.6) 



where, iW = Z^\5 



xe5 ye5= ' •' ' 



(2.7) 



is the self-energy of the droplet 6, which is negative. Moreover, the third sum 



on the r.h.s. of Eq. (2.6) runs over unordered pairs of distinct droplets, and 



W{6,6'):=4J2J2] fp (2.8) 

f ^ f ^ X — V K 
xe(5 ye5' ' ^ ' 

is the droplet-droplet interaction, which is positive. Note that the choice of 

+ boundary conditions implies that all the droplets are closed and within A. 

Our first goal is to get a lower bound on the droplet's self-energy, which is 

suitable for later localization of the energy into small squares of side i, with 



i ^ h*, where h* is the optimal stripe width, see Eq. (2.2). For this purpose, 
given a droplet 6 G T>{A) and the corresponding boundary T{6), we define the 
notion of "bonds facing each other in 6" , in the following way. Let us suppose 
for definiteness that 6 G r(5) is vertical and that it separates a point Xf, G 5 
on its immediate right from a point y^ = x^ — (1, 0) ^ S on its immediate 
left. Consider the bond b' G T{6) such that: (i) b' is vertical; (ii) b' separates 
a point Xfo/ G (5 on its immediate left from a point y^' = Xf,/ + (1,0) ^ 6 on 
its immediate right; (iii) the points x^ and x^/ are at the same height, i.e., 
[x&]2 = [x6']2, and all the points on the same row between them belong to 
S: in other words, x;, -|- (j,0) G 6, for all j = 0, . . . , [s.b']i — [^b]i- We shall 
say that b' faces b in 6, and vice versa. An analogous definition is valid for 
horizontal bonds. Note that in the presence of + boundary conditions all the 
bonds in T{6) come in pairs b,b', facing each other in 6. 

In Appendix p^ we show that the self-energy U{S) can be bounded from 
below as 

m>-y: E E^^^^^^¥f^^^+2^-^^cW+4 E r^^ 

~ ^-^ ^-^ ^-^ n P ^^-^ X — V P 

i=l,2 h&i{&) n^O ' ' {x,y}eP(<5) ' ^ ' 

(2.9) 
where: 

• rj((5) is the subset of r(5) consisting of bonds orthogonal to the i-th 
coordinate direction. 

• db{S) is the distance between b and the bond U facing it in 8. 

• Nc{6) is the number of corners of T{S). 

• V{6) is the set of unordered pairs of distinct sites in 6 such that both 
Cx^y and C^^y cross at least two bonds of T{6). Here C^^y is the path 
on the lattice that goes from x to y consisting of two segments, the 
first horizontal and the second vertical. Similarly, C^^ is the path 



on the lattice that goes from x to y consisting of two segments, the 
first vertical and the second horizontal (note that the two paths can 
coincide, in the case that Xj = yi for some i G {1, 2}). 



The lower bound in Eq. (2.9) is very convenient for localization of the 
energy into small boxes, as shown explicitly in the next section. Let us 
remark that, if desired, the first term on the r.h.s. of this inequality can be 
further bounded from below as 



(2.10) 



3 Localization and minimization 

We introduce a partition of the big box A into squares Q of side £, to be 
optimized in the following. Our purpose is to localize the energy into these 
squares, and to minimize the energy exactly in each small box, thus deriving 
a lower bound on the global energy of the system. Given a droplet configu- 
ration T) and 5 G P, we say that h G r(5) belongs to Q if either it belongs 
to the interior of Q, or it belongs to the boundary of Q and separates a site 
X G 5 n Q from a site y ^ 6. Note that with this definition every bond 
in T{6) belongs to exactly one square Q. The set of bonds b G T{6) be- 
longing to Q will be denoted by Tq{6). The notion that we just introduced 
induces a partition of T{S) into disjoint pieces assigned to different squares: 
T{6) = UqTq{6). Moreover, ii Sq = 6nQ, we define 6^q\ . . . j'^'^^^^'^ to 
be the maximal connected components of 5q, and Fg , . . . , Tg "^ to be the 
portions of Tq{6) belonging to the boundary of ^g , . . . , ^J^*^ , respectively. 
We shall refer to the pair (^g , Fg ) as to a bubble in Q originating from S. 
We shall indicate by Bq{S) the set of bubbles in Q originating from 6, and 
by Bq = Usizx>BQ{6) the total set of bubbles in Q. 

Given /3 = {6,r) G Bq, note that in general F is a union of disjoint 
polygonal curves, each of which can be either closed or open. If one of these 
curves is open, then its endpoints must belong to the boundary of Q. Given 
an endpoint v of an open component of F such that: (1) v is not at a corner 
of Q, (2) the bond 6 G F exiting from v belongs to the boundary of Q; then 
we shall say that F has a "boundary corner" at v. The corners of F belonging 
to the interior of Q will be called "bulk corners" . Moreover, we shall denote 
by Nc{P) the total number of corners of F, i.e., the number of its boundary 



corners plus the number of its bulk corners. Note that 

J2 E NM<N.{S). (3.1) 

This is an inequality (rather than an equality), in general, because S could 
have corners located exactly at the corners of the squares Q. We now derive a 
lower bound on the total energy in terms of a sum of local energies involving 
the bubbles we just introduced. First of all, using Eqs. (2.9) and (3.1), we 
bound the self-energy U{6) from below as 

UiS)>T.{ E UqW) + 1 E W^(A^')}, (3.2) 

Q P&Bq{5) /3,/3'eBQ{<5) 

where the first term on the r.h.s. originates from the first two terms on the 
r.h.s. of (2.9), while the second originates from the last term on the r.h.s. 
of (2.9). The functions Uq and W are defined as follows: if /3 = (5, f), 

^' = ((5',n, 

W^(/3,/3')=4En^' (3.3) 

ye5' 
while 

vm = - E E E """"ijf'" + 2-iiv.OT . (3.4) 

In the last formula, fj is the subset of Fj consisting of bonds orthogonal to 
the i-th coordinate direction, and d^{S) is the distance between h and the 
bond h' e F(5) facing it in 5, if both h and h' belong to f , otherwise it is 
infinite. In a similar manner, we can bound the droplet-droplet interaction 
from below as 

W{5,5')>Y,{ E ^(^"^')}. (3.5) 

/3'gBq(5') 

Inserting Eqs. (3.2)-(3.5) into Eq. ( |2.6[ ) gives 

H^{a^)>J2EQm, (3.6) 

Q 



where 



Eq{Bq)= E [2J|f| + [/Q(^)]+l E ^i^^^')- (3-7) 



2 

p=(lT)eBQ prP'&BQ 



Now consider a bubble (5 = {6, T) such that Nc{f3) > 0, i.e., 6 is not a stripe. 



Proceeding as in Eq. (2.10), we can bound Uq{j3) as 

Uq{P)>-2J,\T\+2'-'2N,{P). (3.8) 

Therefore, 

2J|f I + Uq{^) > r|f I + 2^-^^N,0) . (3.9) 

Note that, in order for T to be very long, the number of corners must be 
sufficiently large: in formulae, 

\T\<2£ + 2iN,{P) . (3.10) 

[The reason is: (a) T (which, in general, is a disjoint union of polygonal 
curves) can have at most two exactly straight lines, and this accounts for the 
2i. (b) Associated with each corner is an ell-shaped open curve, completely 
contained in Q, with the corner at the apex of the curve. The length of this 
curve is at most 2i, and it is clear that the union of all these curves covers 
the remaining part of T. This accounts for the 2iNc{P).] 

If, as we are assuming, NdP) > 0, then NdP) + 1 < 2iVc(/3), so that 

NciP) > ^ . (3.11) 



Inserting this back into Eq. (3.9) gives 



2J|f I + Uq{^) > 2'-'^(l - 4 • 25-Vl^) , (3.12) 

which is positive as soon as \t\£ < 2^^2/4. Therefore, for i shorter than 
2^^2/(4|r|), we can decrease the local energy Eq{Bq) by erasing all the 
bubbles with at least one corner. Denoting by Sq C Bq the subset of Bq 
consisting of bubbles without corners (i.e., consisting of stripes), this means 

that, if £<2i-5/(4|r|), 

Eq{Bq) > Eq{Sq) . (3.13) 

Let now Sq = {f3i, . . . ,(3^}, and let us assume without loss of generality 
(w.l.o.g.) that the stripes (3i = (^j,fj), i = l,...,m, are vertical, and are 
numbered in a way compatible with their order, from left to right. Let us 
also assume w.l.o.g. that Q = [l,i]'^ fl Z^. If m = 1 and /3i = (Q,0), then 
Eq{Sq) = 0. Let us then assume that Fi 7^ 0. Note that the contours Fj 
consist of two vertical parallel lines, for all 2 < i < m — 1. li i = 1, the 
contour Fi can either consist of one or two vertical parallel lines; in the first 

9 



case, 6i = [l,yi] x [1, £] flZ^ for some integer 1 < ii < ^, Fi is the vertical line 
located at the horizontal coordinate yi + ^, and we shall say that Q has — 
boundary conditions on the left; in the second case, 6i = [yo, yi] x [1, £] fi Z^ 
for some integers 1 < yo < yi < i, Ti is the pair of vertical lines located 
at the horizontal coordinates yo — ^,yi + ^, and we shall say that Q has 
+ boundary conditions on the left. Similar definitions are valid for Tm and 
for the boundary conditions on the right. Note that we can always reduce 
ourselves to the case where both the left and right boundary conditions are 
+, up to an error term that is negligible provided that (i ^ h* . In fact, 
suppose that the boundary conditions on the left (say) are — : then we can 
change them to + by erasing the bubble /3i, thus increasing the energy by at 
most \t\L This error term is much smaller than ^^estripes(^) — £^|r|'^P~^^/*^P~'^^ 
if £ ^ /i*. Calling Sq C Sq the set of stripes obtained from Sq after the 
possible erasing of /3i and /3m, we then have 

Eq{Sq)>Eq{Sq) + 2tI. (3.14) 

By construction, Sq consists of k vertical stripes, with k G {m — 2, m — 1, m}, 
whose contours are located at the horizontal coordinates \ < xi < X2 < ■ ■ ■ < 
X2k < ^ + |- We define hi = Xj+i — Xj, with i = 1, . . . ,2k — 1. At this point 
we can utilize the reflection positivity of the kernel |x — y|~^ (see [HI [T7]). 
which leads to the chessboard estimate proved in [121 [131 [IS] • This estimate 
yields the inequality (see Appendix [C] for details) 

2fc-l 

Eq{Sq) > £ ^ /i, [e,{hi) - Chl-H'^] +Ti- Ct-^ , (3.15) 

where es{h) is the specific energy of the periodic striped configuration with 
stripes all of size h, defined in the introduction, and C > is a suitable 



constant. In order to get a lower bound on the r.h.s. of (3.14), we can 
minimize the expression in square brackets over hi < i: 

min [e,{hi) - Ch^Pr'] = estnpes(J)(l + (const.)rVr'/^^"'^)) , (3.16) 
which follows from the explicit expression of es{h) computed in Appendix JA] 



provided that i ':^ h*. Inserting (3.16) into (3.15), and using the fact that 
Zli hi < i, we get 

Eq{Sq) > fesiJ) (1 + (const.)r Vr'/^^"'^) + ri- Ct-^ , (3.17) 

where e^{J) = —Cs{t)\t\^p~'^^^^p~^^ is the optimal striped energy per site in 
the thermodynamic limit. Moreover, the minimum in the first line of the last 



equation is attained at h{i) = h*{l + 0{h*/i)) with h* given by Eq. ( [2.2| ) 

10 



Putting things together, we find that, for i < 2^ 2/(4|r|), 

"^'^^ es(J)(l + 0{r^\r\-'/^P-^^)) . (3.18) 



|A| 

The optimal choice of £ is £ ~ 1^1^^, which gives (recaUing that eo(J) is the 
actual ground state energy per site of our problem): 

^o{J) ^ ^ ^„„„„^ M_|(p-4)/(p-3) 



es(J) 



>1- (const.) |r|^P-4^/^P-^^ . (3.19) 



This proves Eqs. (1.3)-(1.4) and is our final result in two dimensions. In 



three dimensions we can repeat a completely analogous proof, see Appendix 



D} the final result being 

^^ > 1 - (const.) |r|^P-«^/(^P-«^ . (3.20) 



es{J) 

D 

To conclude, let us remark that the proof above also shows that the more 
there are corners, the larger the energy becomes: in formulae, 

H^iaA) - |A|es > c,N, - C^i\A\r''/^''-'^ + \dA\) , (3.21) 

where Nc is the total number of corners associated with a a and ci,Ci > 
are two suitable constants. Therefore, in the ground state, irrespective of the 
boundary conditions, if |A| is large enough, Nc < (const.) |A|r'^/*-'^~^^. In other 
words, by partitioning the macroscopic box into squares of side i' <^ r^'^/i^.-'^)^ 
only a fraction (^t^/(^~^)£'Y of these squares contains a corner of a\, i.e., the 
large majority of these squares are such that the corresponding restriction of 
the ground state is striped or slabbed. A similar argument shows that most 
of these striped/slabbed restrictions consist of stripes or slabs all of a width 
very close to the optimal width h*. 

A Computation of the energy of the optimal 
periodic state 

The specific energy of a periodic striped/slabbed configuration in our two- 
or three-dimensional system is the same as the specific energy of a peri- 
odic striped configuration in an effective one-dimensional system with formal 
Hamiltonian 

H = -jJ2i(^^(^y - 1) + J2^^^^y ~ ^)^(^ ~ ^) ' ^'^•^) 

{x,y) x<y 

11 



where, for all x 7^ 0, 



^(^) = Yl 



nG/ 



■^ (x2 + |n|2W2 • 



(A.2) 



The interaction potential v{x) can be conveniently rewritten as v{x) = V{x) + 
R{x), where 



V{x) 



dy 



-1 (X2 + |y|2)P/2 \x\P-d+^ 7jj,_i (1 + y2)P/2 



dy 



rin- 



Ip-d+l 



\X 



(A.3) 

and R{x) is a rest, which decays to zero at infinity exponentially fast (as 
one can prove by using Poisson's summation formula). The energy of a one- 
dimensional periodic state consisting of blocks all of the same size h and 
alternating sign is straightforward to compute, and the computation gives 
(see m Eq. (17)]): 



e^{h) 



2J 2 r^ ^ . ^ e-" , ah 



r 2 



6 ^ oth 

da^i^{a) ^_^ [l -tanh— ; 



(A.4) 



where yUi,(a) is the inverse Laplace transform of f (x), i.e., the function such 
that v{x) = /g°° da fiv{a)e~°'^ , Wx > 0. Of course, /i^(a) can be rewritten as 
fiv{a) = /iv(Q;) + yUij(Q;), according to the decomposition i;(x) = V{x)+R{x), 
with /iy(a) = YT~^^^d+T)^^~'^ ^^^ /^-r(q^) is zero for a sufficiently small. Plug- 



ging this into Eq. (A.4) and computing the resulting integral asymptotically 
as /i — )■ 00 gives 



es(/i) 



r , A{p) 



fO(- 



1 



;), A4p) 



K„ 



h hP-'^ V/iP-d+2/ ' "V^/ Y{p-d+l) 

Finally, optimizing over h gives 

-{p-d)Ad{p) 



POO 

2P-'^ / daaP~'^-\l-tanha) 
Jo 

(A.5) 



h* = argmines(/i) 
es{J) = e,{h*) = - 



p-d-l 



p — d — 1 



(l + 0(|r|.^=3^)) 



(A.6) 



p — d 



[{p - d)P"iAdip)] 



i/{p-d-i) 



p-d-i(^l + o(|r|p-'«-i)) 



which proves Eqs. (2.1)-(2.2). 
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B Proof of Eq. (K^ 



We start by proving a weaker version of Eq. (2.9), namely 



j=l,2 b^TiiS) nT^O 



mm{\ni\,db{S)} 



nP 



(B.l) 



Later we will show how to improve (B.l) to (2.9). Let us rewrite the droplet 
self-energy as follows: 



U{6) 






nP 



(B.2) 



where A/'n(5) is the number of ways in which n = {rii, 77,2) may occur as the 
difference y — x or x — y with x G 5 and y ^ 6. Let Ti{6) be the subset of 
T{6) orthogonal to the z-th coordinate direction. Our claim is that 



^fn{S)<J2 E min{|n,|,4(5)} 



(B.3) 



from which Eq. (B.l) readily follows. If ni = or ^2 = 0, then the proof of 



(B.3) is elementary, and we leave it to the reader. Let us consider explicitly 



only the case that both rii and ^2 are 7^ 0. We need to define a few geometric 
objects, which are illustrated in Figure |2j 

Consider a pair of points x, y such that x G 5 and y ^ 6. Draw 

-•hv 



the oriented lattice path C^'^l^y that goes from x to y and consists of two 
segments, the first horizontal and the second vertical. Let 61 = 6i(x, y) 
be the first bond in T{6) crossed by C^^^, bi separates a site x;,^ G 6 



from a site y^^ ^ 6. Moreover, let x' 



we define 



be the corner of C^^y] 
or Zb^ (x, y) = x' if bi is between 



Zfej (x, y) = X if bi is between x and x' 
x' and y. This construction allows us to associate the pair (6i,ZbJ with 
(x, y). Similarly, drawing the oriented lattice path C^^y that goes from x 
to y and consists of two segments, the first vertical and the second hori- 
zontal, we can associate with (x, y) a second pair (62,262). By construc- 
tion, in both cases the distance of z^. from y^. is < min{|xj. — yj-\,db-{S)}, 
where ji = 1 if 6j is vertical, and ji = 2 if 6j is horizontal. We write 
•?^(x,y) = {(6i(x,y),Zb^(x,y)(x,y)),(62(x,y),Zfe2(x,y)(x,y))}. Vice versa, if 
we assign an integer vector n 7^ 0, a bond b G Ti{6) separating xj, G 5 from 
y6 ^ 6, and a site z^ G Zb{n, 6) (here Zb{n, S) is the set of allowed locations 
of Zb, namely, is the set of points Zb E 5 belonging to the same column/row 
as b depending on whether b is horizontal/vertical, with the property that 
|zb — yfe| < min{|nj|, db{S)} and all the sites between Zb and x^ belong to 6), 
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Figure 2: An illustration of the geometric objects introduced after Eq. (B.3) 



The grey area is the droplet 6. The two dotted paths connecting x with y 
are C^^^ and C^\y The intersection of the two paths with the boundary 
r((5) defines the two special bonds bi = 6i(x, y) and 62 = &2(x, y). Every 
such bond is associated with a point in 6, denoted by zj,. and located on the 



path C^'^y or C^^y, which can coincide or not with x. 



then the set Qn{b, z^,) = {(x, y) G J-" ^{b, Z5) : x— y G {±n}} has at most two 



elements. This fact immediately implies Eq. (B.3). In fact, if x{condition) 
is the function = 1 when condition is verified, and = otherwise, 



Mr.{5) = 5^x(x-yG{±n}) 

xe(5 

= ^E^(^-y^{±^»E E x((&,z,)G^(x,y)) 



xSiS 

y0 



Zi,e2t,(n,(5) 



\tl E E^((^'y)^^'^(^'^^)) 



(B.4) 






< E E min{|ni|,4((5)} 



where in the last inequality we used the facts that \Qn{b,7.h)\ < 2 and 
\Zb{n,8)\ < min{|ni|,4(5)}. 



Let us now discuss how to improve (B.l) into (2.9). First of all, from its 



proof, it is clear that (B.3) overcounts the pairs in V{5) (for the definition 
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of V{S), see the fourth item after (2.9)). Therefore, we can freely subtract 
from the r.h.s. of (B.3) the additional contribution coming from these pairs, 
so that 



^fniS)<J2 E min{|n,|,4W}-2|P„(5)| 



(B.5) 



where Vn{S) = {{x,y} G V{6) : x - y G {±n}}. Inserting (B.5) into (B.2) 
gives 



i=l,2 beFiiS) n^O 



mm{\ni\,db{6)} 



nP 



E 



1 



{x,y}6P(<5) 



I In ' 



(B.6) 



which is almost what we are after, up to the term in (2.9) proportional to 
Nc{S). In order to get it, let us consider the special case of n such that 
1^1 1 = 1^2! = 1- Note that if |ni| = \n2\ = 1, then Zb{n,6) consists of a single 
point, V6 G T{6). The key remark is that for every bond 6 G r(5) adjacent 
to exactly one corner of T{6), we have 



1 



J2 \On{b,Zb)\<S 



(B.7) 



n: |ni| = |n2|=l 

while for every bond 6 G r(5) adjacent to two corners of T{S) 

1 



Y. \0n{b,Zb)\<2. 



(B.8) 



n: |ni| = |n2|=l 



Of course, in the last two equations z^ is the unique element of Zb{n, 6). Note 
that in general ( |B.7 ) is an inequality (rather than an equality), because the 
corner which b is adjacent to could actually be a "double-corner" like one of 
those in FigjT| rather than a standard one; in fact, if b adjacent to exactly 
one double-corner of T{6), then ^J2n- |ni|=|n2|=i l^n(^' ^^)l = 2. A similar 
comment is valid for Eq. (B.8). 

Using the same rewriting as in Eq. (B.4), together with (B.7)-(B.8), we 



find 



I J2 A/'nW = 7 E E E E^((^'y)^^'^(^'^^)) 



n: |ni|=l, 
|n2|=l 



n: |ni|=l, i=l beri(<5) xG(5 
|n2|=l Zi,G26(n,5) y^S 



^ ^ E E ^i^^i^^i'^^(^)>-^^(^)- (^-9) 



n: |ni|=l, J=l,2 
|n2|=l ber,{5) 
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Moreover, if we also take into account the presence of double- corners, as 



discussed after (B.8), then we can further improve (B.9) into 



\ J2 ^n{S)<l Yl J2'^^i\''^\^d,{6)}-N^{6)- J2 2- 



n: |ni|=l, 
|n2| = l 



n: |ni|=l, i=l,2 

|„2|=i ber,(5) 



{x,y}eP(5): 
|a:i-|/i| = |a;2-?;2|=l 

(B.IO) 



Combining (B.IO) with Eqs. (B.2) and (B.5) finally gives Eq. (2.9). 



C Proof of Eq. (3.15) 



Let Sq = {(3i, . . . ,f3k} be a bubble configuration consisting of k vertical 
stripes, with + boundary conditions on the left and right sides of Q. We 
assume that the bubbles' contours are located at the horizontal coordinates 
^ < xi < X2 < ■ ■ ■ < X2k < ^ + ^, and we let hi = Xj+i — Xj, with i = 
1, . . . , 2k— 1. Given the spin configuration aqinQ = [1, £]^nZ^ corresponding 
to Sq, we can naturally extend it to the strip A2l,£ = [—L + 1, L] x [l,i] flZ^, 
by filling the portions of A2L/ to the left and to the right of Q by + spins; 
we denote the resulting spin configuration by cta^^ ^ . By construction, the 
droplets' boundaries within A2l,£ are still located at xi < ■ • ■ < X2m- 

In terms of these definitions, we can rewrite the energy Eq{Sq) as follows: 



k 



n.P 



l<i<j<k xg5, 



|x-y| 



Ajki~2y y ^^+4 y y-^ 



ye5f\S'i 



|x-y 



\<i<j<k xe^i 



x-y| 



(C.l) 



where in the second line bl = 7? \5i and Si is the infinite vertical strip of 
width /i2i-i containing 5,, i.e.. Si = [{x2i-i, X2i) nZ] x Z. It is convenient to 
rewrite 6i \ Si = Ai U Bi, where Aj = Z^ \ {Si U Aqo/) and Bi = Aqo/ \ Si. 
Correspondingly, we can rewrite: 



EQiSQ) = uki-2 y^ ^^+4 E u^ + E^^(^--i) 

i=l,...,k ' -^ ' l<i<j<k ' -^ ' J=l 

xe(5i, yS-Bi xe5i, ye<5j 

k 

hm i/S(5^A...) + E Mh2.-i) , (C.2) 



L— ^00 



i=l 
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yp 



(C.3) 



where H^^' is the finite volume Hainiltonian (2.4) with periodic bound- 
ary conditions in the horizontal direction and open boundary conditions in 
the vertical direction (of course, the choice of boundary conditions in the 
horizontal direction is arbitrary in the limit L — )■ oo), and 

xe5i ' ' -h2t-^i<xi<o yi,y2>0^ 

y€Ai -i<X2<0 

Note that, for h<(i, 

At{h) = -K + 0{h'^-'P) , (C.4) 

where k is a positive constant independent of d, and h (it coincides with the 
"corner energy" defined in [171 Eq. (3)]). The spin configuration aAa^^ we 
are interested in is quasi-lD, i.e., the value of o'[xi,x2) is independent of X2. 



We shall write a 
Correspondingly, 

where 



(xi,x2) = 0-xi and cTA2jvL,f = ^A2L5 with A2L = [-L + 1, L] n Z. 

^Hr;\aK2.) , (C.5) 



H 



per,0 /~ 



sr'' 



■J 2J i'^xCTx+l 
-L<x<L 



M^ - y) 



1) + X^ (j)i>{x - y){a^ay 

-L<x<y<L 
1 



1) 



-E E 

^ gez m,n=i [{x - y + 2qLf + (m - nf] 



i 



p/2 



(C.6) 



is a one-dimensional spin Hamiltonian with a refiection positive long-range 
interaction and periodic boundary conditions, of the class considered in [121 
[T3] . Therefore, we can apply the chessboard estimate proved e.g. in the 
Appendix of [13]. As a result, using pfS^, Eq- (A5)] and recalling the fact that 
the periodic spin configuration ctaj^ consists of blocks of alternating sign, of 
size hi,. . . , /i2fc-i, h2k, where /i2fc = h2k{L) = 2L + i + xi - X2k, we get 



^:?(-A. 



2k 

> ^ hiee{hi 



(C.7) 



where ee{h) is the energy per site (as computed from H^^' , in the limit 
L — 7- 00) of the infinite periodic configuration consisting of blocks all of the 
same size h, and of alternating sign. Inserting Eqs. (C.4) and (C.7) into 



Eq. (ia2j), we find: 

EQiSg) > 
+ 



2k-l 



ihiCiihi 



K 

2 



O 



hr 



K 

-2 + 



(C.^ 



' lim h2k{L)eAh2k{L)) 
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Now we observe the following: 

ihe^ih) = iheeih) - ^ + O(^) + O(^) , (C.9) 

where es{h) is the specific energy of the infinite periodic striped configuration 
defined in the introduction. Moreover, recalling that \iinL^c>o h2k{L) = +00 



and using (|C.9|) together with (|A.5|), we see that Mnih^oo he^ih) = t + ^ + 

o{e-p). T 



lerefore, for a suitable constant C > 0, 

2k-l 



Eq{Sq) > £ ^ /li [es(/i.) - Chl-n-^] + r£ - C^-^ , (CIO) 



i=l 



which proves Eq. (3.15) 



D Three dimensions 

In this Appendix we adapt the argument spelled out above for two dimen- 
sions to the case of three dimensions, by introducing droplets and contours 
analogous to the two-dimensional ones. Note that now bonds separating a + 
from a — spin are replaced by plaquettes. Droplets now are three-dimensional 
regions whose boundaries are unions of plaquettes. The energy still admits 



the representation (2.6). The first issue to be discussed is the lower bound 



on the self energy of the droplets, which should be replaced by the analogue 



of (2.9), namely 

i=i feer,(<5) riT^o ' ' {x,y}ep(<5) ' -^l 

(D.l) 
where now the label b G Ti{6) is associated with a plaquette of the boundary 
of 5, orthogonal to z-th coordinate direction, and db{6) is the distance between 
b and the plaquette b' facing it in 6. Moreover, Nc{S) is the number of edge 
corners belonging to T{6). By 'edge corner' we mean an edge that is common 
to two orthogonal plaquettes of T{6). Note that an edge corner has length 
1. Finally, V{6) is the set of unordered pairs of distinct sites in 6 such that 
each of the paths C^^y, C^^y and C^^y cross at least two bonds of T{6). 
Here Cl^Xy is the path on the lattice that goes from x to y and consists of 
three segments, the first in coordinate direction i, the second in coordinate 
direction j and the third in coordinate direction k. 



The proof of (D.l) follows the same lines as the proof in Appendix B 



The only relevant differences are the following. When constructing the set 
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J^(x, y) we have to draw the three disjoint lattice paths C^y, C^^y and 
C^lly, so that J-'(x, y) consists of exactly three elements. Similarly, the set 
Gn{b, Zfe) consists of at most three elements. From these considerations, the 



analogue of (B.4) immediately follows. 



The proof of (B.5) is unchanged, and the proof of (B.IO) does not even 



need to be repeated or adapted. Indeed, the analogue of the l.h.s. of (B.IO) 
that we now want to estimate is 



1 

9 


>: 


AfniS) 


1 

+ 9 


^ 


AfniS) 


1 

+ 9 


>; 




ni\ = \n2\ = 

713=0 


=1 




ni|=|?i3| = 
712=0 


= 1 




™3l = |l2| = l 
711=0 



AfniS) 



(D.2) 



Note that the n vectors involved in these sums are all the vectors whose length 
is a/2. The first sum, for example, is really a sum over the contributions 
from horizontal sections of 6, at constant x^; each of these can be estimated 



in exactly the same way as in (B.IO). The same holds for the second and 



third sums above. Putting all these together allows us to estimate (D.2) from 
above by 



1 
2 



Y, J2 ^HH,db{6)}-N,{6)- Y. 2 



(D.3) 



n: |n|=V2 *=1,2,3 



{x,y}eP(<5): 

|x-y|=V2 



which is the desired analogue of (B.IO) 



The next step is localization into boxes of side i. The relevant definitions 
remain unchanged (with certain obvious changes, e.g., the summation over 
z = 1,2 in (3.4) should become i = 1,2,3), and the key estimates (3.7)-(3.9) 
are still valid without alteration. The symbol /3 = {6, T) will still indicate a 
bubble (i.e., a pair consisting of a droplet and its contour; the bars are meant 
to remind the reader that both the droplet and the contour are localized into 
a finite box); similarly, NdP) will still be the total number of corners of F, 
i.e., the number of its boundary corners plus the number of its bulk corners; 



see the lines preceding (3.1 ) 
should be replaced by 



The first estimate to be changed is (3.10), which 



|F| < 2£2 + 2£iV,(/3) 



(D.4) 



The reason is completely analogous to the one explained after (3.10). Insert 



ing (D.4) into (3.9) gives 



2J|f I + UqW) > 2'''^NM{1 - H^ - 2^|r|£) 

If N^^) > 1, then 

2 J|f I + Uq0) > 2'-"^N,0) (1 - 2i+i \r\f) , 



(D.5) 



(D.6) 
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which is positive as soon as i < 2^2^4 |r|^^/^. 

Under this condition, therefore, for the purpose of a lower bound, we can 
erase all the bubbles with one or more corners, and obtain the analogue of 



(3.13). It is at this point that columns are excluded, because a column has 
many edge corners. From this point on the proof is very similar to the one 
of the two-dimensional case: We can assume without loss of generality that 
our bubble configuration of interest consists of a collection of straight slabs. 
Moreover, we may reduce ourselves to + boundary conditions, up to an error 



of the order ri"^, so obtaining the analogue of (3.14), with 2r£^ replacing 2r 



in the right hand side. Now we are in conditions to apply reflection positivity. 



the result being the analogue of (3.15), namely 

2k-l 



i=l 

where now es{h) denotes the energy per site of the periodic slab energy. 
Minimization of this expression under the required constraints on hi and i 
leads to our final result, 

^°^ '' > 1 -(const.)|r|^^ . (D.8) 



^slabs \'^ ) 
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